Cloning transformations in spin networks without external control 
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In this paper we present an approach to quantum cloning with unmodulated spin networks. The 
doner is realized by a proper design of the network and a choice of the coupling between the qubits. 
We show that in the case of phase covariant doner the XY coupling gives the best results. In the 
1 -> 2 cloning we find that the value for the fidelity of the optimal doner is achieved, and values 
comparable to the optimal ones in the general N — > M case can be attained. If a suitable set of 
network symmetries are satisfied, the output fidelity of the clones does not depend on the specific 
choice of the graph. We show that spin network cloning is robust against the presence of static 
imperfections. Moreover, in the presence of noise, it outperforms the conventional approach. In this 
case the fidelity exceeds the corresponding value obtained by quantum gates even for a very small 
amount of noise. Furthermore we show how to use this method to clone qutrits and qudits. By 
means of the Heisenberg coupling it is also possible to implement the universal doner although in 
this case the fidelity is 10% off that of the optimal doner. 

PACS numbers: 03.67.Hk,42.50.-p,03.67.-a 



I. INTRODUCTION 

The no-cloning theorem p| states that it is impossible 
to make perfect copies of an unknown quantum state. At 
variance with the classical world, where it is possible to 
duplicate information faithfully, the unitarity of time evo- 
lution in quantum mechanics does not allow us to build 
a perfect quantum copying machine. This no-go theo- 
rem is at the root of the security of quantum cryptogra- 
phy |2| , since an eavesdropper is unable to copy the infor- 
mation transmitted through a quantum channel without 
disturbing the communication itself. Although perfect 
cloning is not allowed, it is, nevertheless, possible to pro- 
duce several approximate copies of a given state. Several 
works, starting from the seminal paper by Buzek and M. 
Hillery 3] , have been devoted to find the upper bounds to 
the fidelity of approximate cloning transformations com- 
patible with the rules of quantum mechanics. Besides the 
theoretical interest on its own, applications of quantum 
cloning can be found in quantum cryptography, because 
they allow us to derive bounds for the security in quan- 
tum communication Q, in quantum computation, where 
quantum cloning can be used to improve the performance 
of some computational tasks and in the problem of 
state estimation 5]. 

As mentioned above, the efficiency of the cloning trans- 
formations is usually quantified in terms of the fidelity of 
each output cloned state with respect to the input. The 
largest possible fidelity depends on several parameters 
and on the characteristics of the input states. For an 
N — > M doner it depends on the number N of the in- 
put states and on the number M of output copies. It 
also depends on the dimension of the quantum systems 
to be copied. Moreover, the fidelity increases if some 
prior knowledge of the input states is assumed. In the 



universal cloning machine the input state is unknown. A 
better fidelity is achieved, for example, in the phase co- 
variant doner (PCC) where the state is known to lie on 
the equator of the Bloch sphere (in the case of qubits). 
Upper bounds to the fidelity for copying a quantum state 
were obtained in Refs. Q and |(| in the case of universal 
and state dependent cloning respectively. The more gen- 
eral problem of copying N — > M qubits has been also ad- 
dressed 0- The PCC has been proposed in Ref. @- Sev- 
eral protocols for implementing cloning machines have 
been already achieved experimentally 0, fiol ITU Il2| . In 
all the above proposals the cloning device is described 
in terms of quantum gates, or otherwise is based on 
post-selection methods. For example, the quantum net- 
work corresponding to the 1 — > 2 PCC consists of two 
controlled-not (C-NOT) gates together with a controlled 
rotation [l3j . 

The implementation of given tasks by means of quan- 
tum gates is not the the only way to execute the required 
quantum protocols. Recently it has been realized that 
there are situations where it is sufficient to find a proper 
architecture for the qubit network and an appropriate 
form for the coupling between qubits to achieve the de- 
sired task. Under these conditions the execution of a 
quantum protocol is reached by the time evolution of the 
quantum-mechanical system. The only required control 
is on the preparation of the initial state of the network 
and on the read-out after the evolution. This perspec- 
tive is certainly less flexible than the traditional approach 
with quantum gates. Nevertheless it offers great advan- 
tages as it does not require any time modulation for the 
qubits couplings. Moreover, among the reasons for this 
"no-control" approach to quantum information is that 
the system is better isolated from the environment during 
its evolution. This is because there is no active control 
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on the Hamiltonian of the system. Actually, after initial- 
izing the network one needs only to wait for some time 
(to be determined by the particular task) and then read 
the output. Several examples have been provided so far. 
A spin network for quantum computation based only on 
Heisenberg interactions has been proposed [lil ITfij . An- 
other area where this approach is attracting increasing 
attention is quantum communication, where spin chains 
have been proposed as natural candidates of quantum 
channels EE EH EE IS E3, EH - An unknown quan- 
tum state can be prepared at one end of the chain and 
then transferred to the other end by simply employing 
the ability of the chain to propagate the state by means 
of its dynamical evolution. These proposals seem to be 
particularly suited for solid state quantum information, 
where schemes for implementation have already been put 
forward jHlil. 

Stimulated by the above results in quantum commu- 
nication we have studied quantum cloning in this frame- 
work. The main goal is to find a spin network and an 
interaction Hamiltonian such that at the end of its evo- 
lution the initial state of a spin is (imperfectly) copied on 
the state of a suitable set of the remaining spins. In this 
paper we will show that this is indeed possible and we 
will analyze various types of quantum doners based on 
the procedure just described. We will describe a setup 
for the N —> M PCC and we will show that for N = 1 
and M — 2 the spin network cloning (SNC) achieves 
the optimal bound. We will also describe the more gen- 
eral situation of cloning of qudits, i.e. d-level systems. 
An important test is to compare the performance of our 
SNC with the traditional approach using quantum gates. 
We show that in the (unavoidable) presence of noise our 
method is far more robust. Some of the results have been 
already given in Ref. |25|. In the present paper we will 
give many additional details, not contained in Ref. |25| . 
and extend our approach to cloning to several other sit- 
uations. We discuss cloning of qutrits, universal cloning 
machines, and optimization of the model Hamiltonian 
just to mention few extentions. 

The paper is organized as follows. In Section [H] we 
review the basic properties of approximate cloning show- 
ing the theoretical optimal bounds. In Section IIIII wc 
present the models and the networks topologies consid- 
ered in this work. In Sections II VI and lv| wc briefly review 
and extend the results obtained in |25j. These sections 
concern the spin network model to implement the 1-tM 
and N — » M phase covariant cloning transformations re- 
spectively. In addition to a more detailed discussion, as 
compared to [25j], here we present a detailed analysis of 
the role of static imperfections. We also optimize the 
cloning protocol over the space of a large class of model 
Hamiltonians which includes the XY and Heisenberg as 
limiting cases. 

The effects of noise, included in a fully quantum me- 
chanical approach, are analyzed in Section IVII where we 
compare our cloning setup with cloning machines based 
on a gate design. In Sect. IVIII we study the possibility 



of achieving universal cloning with the spin network ap- 
proach. In Section rVllll we generalize the SNC for qutrits 
and qudits. Finally, in Section Hxl we propose a simple 
Josephson junctions circuit that realizes the protocol and 
in^Jwc summarize the main results and present our con- 
clusions. 



II. OPTIMAL FIDELITIES FOR QUANTUM 
CLONING 

Most of this paper deals with the case of PCC. We will 
therefore devote this section to a brief summary of the 
results known so far for the optimal fidelity achievable in 
this case. 

We start our discussion by considering quantum 
cloning of qubits, whose Hubert space is spanned by the 
basis states |0) and |1). The most general state of a 
qubit can be parametrized by the angles tp) on the 
Bloch sphere as follows 

|^=co S ^|0)+e^sm||l). (1) 

Quantum cloning was first analyzed Q , where the 1 — > 2 
universal quantum cloning machine (UQCM) was intro- 
duced. We remind that the fidelity of a UQCM does not 
depend on (•&, <p) i.e. it is the same for all possible input 
states. As already mentioned, the quality of the doner 
is quantified by means of the fidelity T of each output 
copy, described by the density operator p, with respect 
to the original state \ip) 

^ =<##>■ ( 2 ) 

The value of the optimal fidelity is achieved by maxi- 
mizing T over all possible cloning transformations. The 
result for the 1 -> 2 UQCM is T = 5/6 ~ 0.83 0, @|. 
The general form of the optimal transformation, which 
requires an auxiliary qubit, has been explicitly obtained 
in Ref. @. 

When the initial state is known to be in a given sub- 
set of the Bloch sphere, the value of the optimal fidelity 
generally increases. For example, in Ref. Q cloning of 
just two non orthogonal states is studied and it is shown 
that the fidelity in this case is greater than that for the 
UQCM. The reason is that now some prior knowledge 
information on the input state is available. Another im- 
portant class of transformations, which will be largely 
analyzed in the present paper, is the so called phase co- 
variant cloning. In this type of doner the fidelity is opti- 
mised equally well for all states belonging to the equator 
of the Bloch sphere: 

ty) = -L(\0) + e i *\l}) , (3) 

where <p e [0,27r]. The optimal transformation for the 
1 — > 2 PCC was found in yj. The corresponding fidelity 
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is given by 



T = ~ + -4= ^ 0.854 
2 V8 



(4) 



In the N —> M case of PCC, the optimal fidelities 
were derived in Ref. [2(|. For the 1 — > M PCC case they 
read H| 
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Optimal cloning has also been studied in higher di- 
mensions. The universal case was discussed in Refs. 1291 . 
Cloning of qutrits was specifically treated in Refs. [HI l30j| . 
where the optimal fidelity for cloning some classes of 
states was derived. In particular the double-phase co- 
variant symmetric 1 — > 2 doner , which is optimized for 
input states of the form 
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where (p\ and ip2 are independent phases and the states 
|0) , 1 1) , |2) form a basis for a qutrit, was analyzed 5] . The 
optimal 1 — > M fidelity in the case of M = 3k + 1 (with 
k positive integer) is given by the simple expression |26| : 
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FIG. 1: a. Circuit implementing 1 — ► 2 PCC for qubits and 
experimentally realized in l33il . b. Circuit implementing 1 — > 

3 PCC 0. 



This is not the only way to perform quantum cloning 
and quantum computation protocols in general. In the 
following sections we describe how to accomplish the de- 
sired task using the evolution of a spin network. The ex- 
pressions for the fidelities given above will be compared 
to the ones derived by means of the approach proposed 
in Ref. |25j and in the present paper. 



III. THE SPIN NETWORK CLONING 
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For qudits, namely quantum systems with finite di- 
mension d, the PCC is optimized for states of the form 



d-i 



where ipo = 0, while <£>i,y>2, 
phases. The fidelity for 1 — > 
given by [3l| 



10 



(9) 
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. . , fd-i are independent 
2 PCC in dimension d is 
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More recently the general N — > M case of PCC for 
qudits was analysed |32j, where explicit simple solutions 
were obtained for a number of output copies given by 
M = kd + N, with k positive integer. 

The quantum transformations corresponding to the 
cloning machines described above can be implemented by 
suitably designed quantum circuits. For example, in the 
1 — »• 2 and 1 — * 3 cases the circuits implementing PCC 
were derived in Refs. and respectively and they 
are shown in FigE] (note that in these cases no auxiliary 



qubits are needed) . In Fig|TJ>. we defined R 
and #1 = $3 = ir/8 and $2 
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In this section we show how quantum cloning can be 
implemented using a spin network. First let us discuss 
our Hamiltonian model. We start with a fairly general 
model defined as 



where a l x y z are the Pauli matrices corresponding to the 
j-th site, Jij are the exchange couplings defined on the 
links joining the sites i and j and B is an externally ap- 
plied magnetic field. The anisotropy parameter A ranges 
from (XY Model) to 1 (Heisenberg Model). This model 
is named as the XX Z Model. We discuss separately the 
two limiting cases A = 1 and A = 0. It turns out that for 
PCC the A = case leads always to the highest fidelity. 
Given the model of Eq. the fidelity is maximized 

over B/J and Jt. We defined i?( M ) and t^ M ^ the values 
of the parameters leading to the optimal solution. Notice 
that the total angular momentum as well as its z compo- 
nent are always constants of motion independently of the 
topology of the network. In all the cases we consider in 
this work, the couplings Jy ^ only for nearest neigh- 
bors sites To specify which couplings are non-zero 
one has to define the topology of the spin network. 

For the 1 — > M PCC we choose M + 1 spins in a star 
configuration (see fig. The central spin labeled by 

1 is initialized in the input state while the remaining M 
spins are the blank qubits and are initialized to the state 
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|0) if < i9 < tt/2 and |1) if tt/2 < ■d < n. For this 
network — J only if one of the two linked sites is the 
central one. This configuration has also been studied in 
a different context [3^. For the 1 — *■ M case we have 
considered also other types of networks. These are rep- 
resented in FigEt- The original state is placed on the 
top of the tree while the blank qubits are on the lowest 
level. The intermediary spins are ancillae. Each graph 
is characterized by the number k of links departing from 
each site and the number j of intermediate levels between 
the top and the blank qubits level. The number of blank 
qubits can be obtained from j and k as M — fc J+1 . Notice 
that for this class of graphs a symmetry property holds: 
the global state of the blank qubits is invariant under any 
permutation of the M qubits. For the N — ► M PCC we 
have considered a generalization of the star network (see 
FigEt). It consists of a star with N centers and M tips 
so no auxiliary qubits are present. Also this network is 
permutation invariant. 

The cases defined above are not the most general net- 
works and/or model Hamiltonian conceivable. Since the 
fidelity must be maximized over the parameters of the 
Hamiltonian as well as over the network topology one 
may wonder whether it is sufficient to consider only the 
cases introduced above. In Section ITVB I we partially an- 
swer to this question by considering the more general 
configuration for the 1 — > 3 case containing 4 spins and 
we believe to have found the best possible scenario for 
the SNC. As far as the choice of the model Hamilto- 
nian is concerned, the symmetry in the XY~-plane is sug- 
gested by the phase covariance requirement for a PCC. 
We checked that the Ising model in transverse field, which 
is not phase covariant, gives poorer results for the fi- 
delity. In principle one should also explore the possibility 
of multi-bit couplings, but we did not considered this (in 
principle interesting) situation. Multi-bit couplings are 
much more difficult to achieve experimentally and at the 
end of this work we want to propose to implement our 
scheme using Josephson nanocircuits, where two-qubit 
couplings with XY-symmetry are easy to realize. 



IV. 



M PCC CLONING 



In this case the Hamiltonian can be easily diag- 
onalized because it can be mapped to the problem of a 
spin- 1/2 interacting with a spin-M/2 



H\ — J (SlS x 



Sy Sy 



\SlS z )+B{Sl + S z ) . (12) 



where we have defined S 1 = cr 1 /2 and S = 1/2 ^ 2 a 1 . 
The operators S obey the usual commutation relations 
for an angular momentum operator. Notice that the 
modulus and the z component of the total angular mo- 
mentum F = S 1 + S commute with the Hamiltonian. 
Thus the evolution is invariant under rotations around 
the z axis. This property automatically makes our model 
a PCC. 
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FIG. 2: Different topologies for N — > M doner: a) Spin star 
network for 1 — > M doner, b) Generic graph for the 1 — > M 
doner with j intermediate steps and k links departing from 
each vertex, c) Spin network for the N — > M cloning 



For < i) < tt/2 (the other case is equivalent) the 
initial state of the star is: 



|*(0))=a|0)+j8|l) 



(13) 



where a = cos | and j3 = e l(p sin ^ . Here we use the 
convention that |0) = |0^_0) and \j) = | ■ ■ ■ 1 ■ ■ ■ ) 

M+l M+l 

where only the site jth is in the state |1). 

Because of the conservation of the total angular mo- 
mentum the state of the star at time t will be a linear 
combination of |0) and \j) i.e. states with only one qubit 
in state |1). The state at time t can thus be written in 
the form (apart from a global phase factor): 



M+l 



(14) 



J=2 



where the coefficients (3\ (t) and 02 if) depend on the par- 
ticular choice of the Hamiltonian. In order to calculate 
the fidelity of the clones we need the expression for the 
reduced density matrix of one site. For symmetry reasons 
this is independent on the site chosen, the result being 
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The fidelity for the SNC is 
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where the coefficients (t) depend explicitly on the cho- 
sen model (Heisenberg or XY in this case). 



Heisenberg model - Let us start with the Heisenberg 
model (A = 1) and let B — (we checked that a finite 
external magnetic field is not necessary to achieve the 
maximum fidelity) . The Hamiltonian can be rewritten in 

the form H — JS 1 ■ S and using S 1 ■ S = 1/2(F 2 - 
S 2 — Sf) one finds that the eigenenergies are given 
by E(F, S, S 1 ) = i [F(F + 1) - S(S + 1) - S^S 1 + 1)] 
where F, S, S 1 are the quantum numbers associated to 
the corresponding operators. The eigenvectors can be 
found in terms of the Clebsch-Gordan coefficients. The 
results for the coefficients /3j(t) are: 



ft (*) = 
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(17) 
(18) 



where S = M/2. 

The maximum value for the fidelity 



ft = 



4 + M (3 + M) + (M - 1) [(3 + M) cos <d - cos 2i?] 

2(1 + Mf 

(19) 

is obtained for the parameters Jv> M > = 2tt/(M + 1). 

XY model - Now let us turn our attention to the XY 
model. Solving the eigenvalue problem as in |35| one 
finds 
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The fidelity is maximized when B^ M ' / J 
JtM = tt/VM: 
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(22) 

M/2 and 



(1 + VM) 2 - (2 - 2Af)cos?9+ (l - VM) cos2i? 
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(23) 

For (A = 0) the fidelity is always greater than for (A = 
1). Let us analyze the previous results in two important 
cases. First let us discuss the case M = 2 for arbitrary ■& 
(see FigEJ): 



Ti = 



14 + 5cosi9 - cos 2$ 
18 



(24) 



F Q = i(5 + 2costf + cos2tf + 2\/2sin 2 tf) (25) 
8 

The fidelity Tq coincides with the fidelity for the 1 — > 2 
PCC 36] i.e. the SNC saturates the optimal bound for 
the 1 — > 2 PCC. Second let us consider d — tt/2 and 
arbitrary M (see FigQJ: 
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FIG. 3: The fidelity T\ for M = 2 versus 1? for the XY (solid) 
and Heisenberg (dashed) model are shown. Notice that the 
opti 




FIG. 4: The fidelity T for PCC (circle), XY (diamond) and 
Heisenberg (triangle) as functions of M for $ = tt/2. 



For M > 2, J^x is always smaller than the optimal fi- 
delity given in Ref.|2^|. Also in this case the XY model 
is better suited for quantum cloning as compared to the 
Heisenberg case. Although for generic M the SNC does 
not saturate the optimal bound, there is a very appealing 
feature of this methods which makes it interesting also in 
this case. The time required to clone the state decreases 
with M. This implies that, in the presence of noise, SNC 
may be competitive with the quantum circuit approach, 
where the number of gates are expected to increase with 
M. We analyze this point in Section IVTl 

Recently a PCC with the star configuration has been 
proposed also for a multi-qubit cavity j^]. In this pro- 
posal the central spin is replaced by a bosonic mode of 
the cavity. By restricting the dynamics in the subspace 
with only one excitation (one excited qubit or one photon 
in the cavity) the Hamiltonian is equivalent to the XY 
spin star network considered here. Indeed the optimal 
fidelities coincide with Tq, Ea. H27fl . 

All the results discussed so far have been obtained for 
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the star network. Obviously this is not the only choice 
which fulfills the symmetries of a quantum cloning net- 
work. In general one should also consider more general 
topologies and understand to what extent the fidelity 
depends on the topology. We analyzed this issue by 
studying the fidelity for the XY model and i9 — ir/2 for 
M < 32 for the graph b of FigEl (the fidelity for Heisen- 
berg model in this case is much worse than in the star 
configuration). We conclude that the maximum fidelity 
obtained does not depend on the chosen graph. 



A. Imperfections 

To assess the robustness of our protocol, it is impor- 
tant to analyze the effect of static imperfections in the 
network. In a nanofabricated network, as for example 
with Joscphson nanocircuits, one may expect small varia- 
tions in the qubit couplings. Here we analyze the 1 — > M 
cloning assuming that the couplings have a certain 
degree of randomness. For each configuration of disorder 
Ju are assigned in an interval of amplitude 2e centered 
around J = 1 with a uniform distribution. First we study 
the case of uncorrelated disorder in different links. The 
values of B and t are chosen to be the optimal values 
of the ideal situation. For a given configuration of the 
couplings the fidelities of each of the clones are different 
due to the different coupling with the central spin. Only 
the average fidelity is again symmetric under permuta- 
tion among the clones. We averaged the fidelity over the 
M sites and over 10 3 realization of disorder. For e = I CP 1 
and M < 10 the mean fidelity decreases by just less than 
0.2% of the optimal value. It is important to stress that 
the effect of imperfections is quite weak on the average 
fidelity. This is because for certain values of Jy, even if 
the fidelity of a particular site can become much larger 
than the fidelity in the absence of disorder, at the same 
time for the same parameters the fidelity in other sites 
is very small and the average fidelity is weakly affected 
by imperfections. In figure [S] we show the fidelity for the 
1 — > 2 SNC with imperfections as a function of the toler- 
ance e. We study also the case with correlations between 
the signs of nearest neighbor bonds: the probability of 
equal signs (Ju — J)(Ju+i — J) > is proportional to 
fx G [—1; 1]. The uncorrelated results are recovered for 
fi = 0. As expected this type of disorder is more destruc- 
tive as shown in figure [SJ 



B. Optimal network Hamiltonian for 1 — > 3 PCC 

As shown in FigH the 1 -> 2 SNC saturates the PCC 
optimal bound. However this is not the case for M > 2, 
at least for the network topologies considered up to now. 
One may wonder whether a different choice for the net- 
work could allow to approach the optimal fidelity. In 
order to understand this point we studied the simplest 
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FIG. 5: Mean fidelity for the 1 — > 2 case with static imper- 
fections as a function of the tolerance e with /i = (filled 
squares) and = 0.5 (empty circles). 



non trivial case namely M = 3 and considered the tetra- 
hedron network shown in FigO 




FIG. 6: The tetrahedron network analyzed for the 1 — > 3 
cloning. 

We concentrated on the general anisotropic XX Z 
model presented in 1|11[) in which the local magnetic field 
and the couplings between the central spin and the blank 
spin can be different. This is the most general Hamilto- 
nian for 4 spins that fulfills the symmetry and covariance 
property. For this general model we maximized analyt- 
ically the on-site fidelity diagonalizing the correspond- 
ing Hamiltonian. We found that the maximum fidelity 
exactly coincides with that found with the simple star 
configuration. It is thus demonstrated that, at least for 
M < 3, the star configuration is the optimal network for 
cloning. 

It is however important to stress that, given the trans- 
formation for the optimal PCC, it is always possible to 
find a Hamiltonian that generates this transformation 
during the dynamical evolution. Therefore, at least in 
principle, one should be able to saturate the optimal 
value by including other terms in the Hamiltonian (multi- 
spin coupling for example) . On purpose we chose to limit 
ourselves to a fairly general model which however can be 
realized experimentally. 
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M 


Tpcc 




Jt c 


J/B 


Jt c (KT 2 ) 


2 
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0.941 


0.938 


1516.0 


39.5 


2.9 
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0.933 


0.889 


53.1 


4.5 


5.3 
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0.912 


0.853 


774.1 


12.8 


2.7 
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0.908 


0.825 


563.4 


28.7 


2.8 
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0.898 


0.804 


156.0 


40 


4.9 


2 


8 


0.895 


0.786 


116.6 


29.9 


6.9 


3 


4 


0.973 


0.967 


2201.6 


47.5 


111.8 


3 


5 


0.970 


0.931 


1585.5 


33.1 


19.6 


3 


6 


0.956 


0.905 


8.3 


10.6 


8.3 


3 


7 


0.954 


0.875 


8.1 


3.4 


7.9 



TABLE I: The maximum fidelity T for N -> M for the net- 
work of Fig|2J:. Tpcc is the optimal fidelity for the PCC [2r|. 
Column 5 (6) reports the corresponding evolution time t c (in- 
teraction strength J). Column 7 reports the time t c (e = 10 -2 ) 
at which the fidelity reaches the value T a bs ~ 10 -2 . The re- 
sults refer to the XY model (A = 0). The value T is found 
by numerical maximization in the intervals B/J £ [0.01; 10] 
for N + M < 10 and Jt G [0; 5 • 10 3 ]. 



V. N -> M PCC CLONING 

In this section we discuss the generalization of the 
SNC to the N > 1 case. The suitable network to ac- 
complish this task is depicted in Fig|5J:. The model can 
be mapped to the problem of the interaction between 
two higher dimensional spins, N/2 and M/2 respectively. 
Since we did not succeed in finding the analytic solution 
to the problem (for example for 2 — > 8 the relevant sub- 
space has dimension 56), we simulated it numerically. 
We have simulated the evolution of the network in the 
range B/J e [0.01; 10] and tJ<5- 10 3 . We found the 
absolute maximum of the fidelity T a bs in this interval. 
The result of this maximization is summarized in Table 
d for several values of N and M. We also calculated the 
time to reach a value of fidelity slightly lower than J- a bs- 
The time needed to reach T a bs — S, 5 <C 1, is greatly re- 
duced. Indeed the fidelity is a quasi periodic function of 
time approaching several times values very close to T a bs ■ 
In Table [I] both the absolute maximum T a bs (column 4) 
in the chosen interval and the time t c (5 = 10 -2 ) (last 
column in the table) are shown. 



VI. QUANTUM CLONING IN THE PRESENCE 
OF NOISE 



So far we have described the unitary evolution of iso- 
lated spin networks. Real systems however are always 
coupled to an environment which destroys their coher- 
ence. In this section we will try to understand the effect 
of noise on the SNC. We will also compare the perfor- 
mances of quantum cloning machines implemented with 
spin networks and with quantum circuits using the same 
Hamiltonian. The effect of the environment can be mod- 



eled in different ways. One is to add classical fluctuations 
to the external magnetic field B or the coupling J. These 
random fluctuations can be either time independent or 
stationary stochastic processes. In both cases one can 
define an effective field variance A and average the re- 
sulting fidelity. In Fig. 0we compare the fidelity T\-^i 
and d = tt/2 as a function of A for the XY"-model with 
the optimal average values for fluctuating J (solid) and 
B (dashed). The probability distributions are chosen to 
be Gaussian. Note that the fidelity is more sensitive to 
fluctuations of B. 




FIG. 7: Fidelity for equatorial qubits in the XY model with 
a classical fluctuating field. T is plotted as a function of the 
variance A for fluctuating J (solid) and B (dashed). 



However there are situations in which the environment 
cannot be modeled as classical noise and one has to use 
a fully quantum mechanical description. Following the 
standard approach, we model the effects of a quantum 
environment by coupling the spin network to a bosonic 
bath. Then we describe the time evolution for the re- 
duced density matrix of the spin system alone, after trac- 
ing out the bath degrees of freedom in terms of a master 
equation [37° 



The Hamiltonian for the whole system is 

(28) 
(29) 



H = Hs + Hr + Hj 

M+l 

Hi = Y,ll X ^ k )<[ a \W + < k ) 

i=l k 
M+l 

H R = I>*(fc) ot(fc)oi(fc) 

i=l k 



(30) 



where H$ is the spin Hamiltonian defined in Eq. 1)11(1 . 
The model is presented for generic M but we will discuss 
the results only for M = 2 and M = 3. We suppose 
that each spin is coupled to a different bath, labeled by 
i, and that all baths are independent, u>i(k) and Aj(fc) 
are the frequency and the coupling constant of the fcth 
mode of the ith bath. It is convenient to define the op- 
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the environment 



erator E t = £) fc Ai(fe) a\{k) + ai{k) 
operator to which the system is coupled. 

The master equation in the basis of eigenstates of Hs 
can be written as: 



dt 



Pab 



^ Tt-abcd Red 



(31) 



'ibcd 



where the indexes a, b, c, d run over the energy eigenstates 
and IZabcd is the so called Bloch-Redfield tensor in the 
interaction picture: 

/>oo 

TZ abcd =J2 dr {G i (r)S> cd + G,(-r)E< cd ) (32) 
„■ Jo 



where 



Id&e" 



and 



(33) 



Kbcd = *ac£(4M4W^ - H) a M) db e^ 
n 

(34) 

with (<J l z )ab = (a\al\b). The function G(t) is the correla- 
tion function of the environment operators in the inter- 
action picture: 



G 1 {t) = Tr PF Ei(T)Ei(p) 



(35) 



The functions G^(r) can be related to the spectral density 
of the bath through 



[Giir)}^ = 2JV 4 (w)J i ( Cl ;) 



(36) 



where [■] indicates the Fourier transform. In Eg. 136(1 
Ni(u>) = (e 13 ^ — l) -1 is the mean occupation number 
of the u>i mode at temperature T — and Ji(u>) — 
7r^ t |Ai(wfc)| S(uj — uJk) is the spectral density. We 
suppose that the bath is Ohmic, as often encountered 
in several situations, i.e. J(uj) has a simple linear depen- 
dency at low frequencies up to some cut-off: 



Ji(u) 



-auie 



j/u>c 



(37) 



The parameter a represents the strength of the noise and 
luc is the cut-off frequency. 

In order to to compare SNC with traditional quantum 
cloning machines we have to consider a specific system 
where the required gates are performed. Obviously this 
can be done in several different ways: we choose the XY 
Hamiltonian as the model system for both schemes. In 
particular we compare the two methods for M = 2 and 
M = 3 equatorial qubits. For the quantum circuit ap- 
proach quantum gates are implemented by a time depen- 
dent Hamiltonian. It has been shown |3a H3 that the 
XY Hamiltonian is sufficient to implement both one and 
two-qubit gates. The elementary two-qubit gate is the 
iSWAP: 



U(iSWAP) 



( 1 \ 

i 

1 

1/ 



(38) 



It can be obtained turning on an XY interaction between 
the two qubits without external magnetic field and letting 
them interact for Jt = 7r/4. By applying the iSWAP gate 
twice, the CNOT operation can be constructed 



iSWAP 



iSWAP 



MS) 



Mi) 



MS) 



FIG. 8: Circuit implementing the CNOT from the iSWAP. 
This circuit is used to implement the quantum cloning by 
means of gates. 



This means that we need two two-qubit operations for 
each CNOT. We simulated the circuits shown in Fig. Q] 
for M = 2 and M — 3 in the presence of noise and 
we calculated the corresponding fidelities. We neglected 
the effect of noise during single qubit operations. This 
is equivalent to assume that the time needed to perform 



this gates is much smaller than the typical decoherence 
time. The results are shown in Fig|5] and Fig^| The 
fidelity for the quantum gates (squares) and that for the 
SNC (circles) are compared as functions of the coupling 
parameter a. Even for small a the fidelity for the cir- 
cuits is much worse than that for the network. Notice 
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0.86 




FIG. 9: 1 — > 2 cloning. Comparison of the fidelity To obtained 
by the spin network method and the quantum circuit [XY 
interaction) discussed in [33l . l34|l in the presence of an external 
quantum noise. Circles and squares refer to the network and 
gates case respectively [•& = 7r/2.). The parameters for the 
environment are /3 = 10/ J and loq = 10 4 J. 




'• 5 0.01 0.02 0.03 0.04 0.05 



a 

FIG. 10: The same as in FigEJfor the 1 -> 3 case. 



environment. 



VII. THE UNIVERSAL CLONER WITH SPIN 
NETWORKS 

It would be desirable to implement also a universal 
quantum doner by the same method illustrated here. In 
this section we briefly report our attempt to implement 
the 1 -> 2 universal doner. In the previous sections we 
demonstrated that for the models presented the fidelity 
is invariant on ip (phase covariance) but still depends on 
i9. This axial symmetry relies on the selection of the 
z-axis for the initialization of the blank spins. In or- 
der to perform a universal doner we need a spherical 
symmetry. This means that both the Hamiltonian and 
the initial state must be isotropic. The first condition is 
fulfilled using the Heisenberg interaction without static 
magnetic field that would break the spherical symmetry. 
The second requirement can be obtained using for the ini- 
tial state of the blank qubits a completely random state. 
In other words the complete state of the network (initial 
state + blanks) is 

p(0) = <8> -1 • 

The maximum fidelity is obtained for Jt = ^ and has 
the value 

T = 13/18 ~ 0.72 

that has to be compared with the value 5/6 ~ 0.83 of 
the optimal universal doner Q. Our model is the most 
general time independent network containing three spins 
and fulfills the required conditions. 



VIII. QUANTUM CLONING OF QUTRITS AND 
QUDITS 



that for M = 3, though without noise [a = 0) the SNC 
fidelity is lower than the ideal one, for a > a* = 2.5 • 10~ 3 
the situation is reversed. This shows that our scheme 
is more efficient than the one based on quantum gates. 
Moreover for M > 3 the time required for quantum cir- 
cuit PCC grows with increasing M while, as discussed 
previously, the optimal t^ M > of the SNC decreases with 
M. This suggests that our proposal is even more efficient 
for growing M. Changing the model does not affect these 
results. Indeed the time required to perform a CNOT us- 
ing Heisenberg or Ising interactions is just half the time 
required for the XY model. 

We also believe that in a real implementation the effect 
of noise on our system can be very small compared to the 
that acting on a quantum circuit. This is because during 
the evolution the spin network can be isolated from the 



Spin network cloning technique can be generalized to 
qutrits and qudits. This is what we discuss in this Section 
starting, for simplicity, with the qutrit case. The cloning 
of qudits is a straightforward generalization. Our task is 
to find an interaction Hamiltonian between qutrits able 
to generate a time evolution as close as possible to the 
cloning transformation. One obvious generalization of 
the qubit case is to consider qutrits as spin-1 systems. 
In this picture the three basis states could be the eigen- 
states of the angular momentum with z component (- 
1,0,1). The natural interaction Hamiltonian would then 
be the Heisenberg or the XY interaction 

H I = J lj S i -Si or Hi = JijiS^S^ + S]rS Y ) (39) 

Alternatively one can think to use the state of physical 
qubits to encode the qutrits. Such an encoding, originally 
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proposed in a different context jljj, uses three qubits to 
encode one single logical qutrit: 

|0) L = |001) 
|1>l - |010) 
\2)l = |100) 

In Ref . 38| it is shown that this encoding, together with 
a time-dependent XY interaction, is universal for quan- 
tum computing with qutrits. In our work however we 
have restricted ourselves to the use of time-independent 
interactions with a suitable design of the spin network. 
For the qubit case the XY interaction is able to swap 
two spins. We know that this is the key to clone qubits 
and so one could try a similar approach also for qutrits. 
However, for higher spin, Hund's rule forbids the swap- 
ping. For this reason we have turned our attention to the 
encoded qutrits to see if swapping is possible. It is simple 
to show that the network depicted in Fiell II satisfies our 
requirements. In the arrangement each dot represents a 

Qutrit 2 Qutrit 1 Qutrit 3 




FIG. 11: Arrangement of the network for qutrits. Each dot 
represents a spin and an ellipse encloses each logical qutrit. A 
line connecting two dots means that the corresponding spins 
interact via XY model. 

spin and three dots inside an ellipse correspond to an en- 
coded qutrit. A static magnetic held A pointing in the 
z direction is applied to the hrst spin. A line connecting 
two dots means that they interact via an XY interaction 
with amplitude J. It can be easily checked that for a sin- 
gle couple of qutrits the exchange processes are possible. 
This network is the generalization of the spin star that 
we analyzed before in which a single qutrit interacts with 
the others. It is easily generalized for the 1 — ► M case 
using three spin stars. The single qutrit Hamiltonian is 
realized applying magnetic holds to the physical qubits. 

In analogy with the qubit doner we will prepare qutrit 
1 in the original state 



= a|0) L +/3|l)L+7|2)i 



(40) 



and initialize the other qutrits in a blank state, for ex- 
ample |0)i. Now due to the interactions the state will 



evolve in a restricted subspace of the Hilbert space: 

Mt)> = a|000) L + /3 1 |100) L + /3 2 |010) L + /33|001) L 
+ 7 i|200)l+72|020)l+73|002) l (41) 

To hnd the fidelity of the clones with respect to the state 
of Ea. H40fl we need the reduced density matrix of one of 
the clones (for example the third). The result, in the 
basis (|0) l ,|1)l,|2) l ), is 



P3 



(l 



\ 



i&r- 

a*f3 3 
a*7 3 



l7 3 | 2 aP* 



<*7 3 

&73 
/?3 73 |T3 r 



(42) 



In order to find the coefficients (3i(t) and Ji(t) we have 
to diagonalize the Hamiltonian. We consider the double 
PCC of Eq. J7J): our model is automatically invariant on 
<fi because there is no preferred direction in the space of 
the qutrits. The maximum fidelity achievable with SNC 
is: 



2^2 



0.759 



(43) 



This value has been obtained with A/ J = \j\f2 and 
Jt = 7T/V2- Note that this value is very close to the 
optimal one and the difference is only 2 • 10~ 3 . 




FIG. 12: The optimal (square) and the SNC (diamond) fideli- 
ties for 1 — >• 2 PCC in d dimensions are compared. 



We calculated also the fidelity for the 1 — ► M cloning 
of qutrits using the star configuration. The maximum 
fidelity is: 



T = 



2 + 4VM + 3M 
9M 



(44) 



obtained for the same value of the star configuration of 
qubits (J^ M ) = tt/VAT and /J = VM/2). 

The generalization to qudits is straightforward. Fol- 
lowing the same approach we encode qudits using d 
qubits to encode each qudit. After some algebra one finds 
the general expression for the PCC in d dimensions. The 
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values f( M ) and £?( M ) are independent from d and the 
expression for the fidelity is: 



>2.d 



(d- l)(d + 2V2) 
2d 2 



(45) 



In Fig^l the optimal and SNC fidelities are compared. 
As we can see, the fidelity of the spin network implemen- 
tation is very close to the ideal one. 



IX. IMPLEMENTATION WITH JOSEPHSON 
NANOCIRCUITS 

The final section of this work is devoted to the possi- 
bility of implementing spin network cloning in solid-state 
devices. Besides the great interest in solid state quantum 
information, nanofabricated devices offer great flexibility 
in the design and allow to realize the graphs represented 
in Fig|2| We analyze the implementation with Josephson 
nanocircuits which are currently considered among the 
most promising candidates as building blocks of quantum 
information processors |4ll 14^ . Here we discuss only the 
1 — * 2 cloning for qubits. The generalization to the other 
cases is straightforward. 



a) 



b) 



FIG. 13: a) A sketch of the charge qubit. It consists of a 
superconducting electron box formed with an applied gate 
voltage V x - The device operates in the charging regime, i.e. 
the Josephson couplings Ej of the junction (crossed box in 
the figure) is much smaller than the charging energy, b) Im- 
plementation of the 1 — * 2 spin network cloning by means of 
Josephson qubits. The unknown state to be cloned in stored 
in the central qubit c while the blank qubits u and d are the 
ones where the state is cloned. The coupling between the 
qubits is via the Josephson junctions of coupling energy Jk- 



In the charge regime a J osep hson qubit can be realized 
using a Cooper pair box |4l| (see FiglT3h). the logical 
state is characterized by the box having zero or one excess 
charge. Among the various ways to couple charge qubits, 
in order to implement SNC the qubits should be coupled 
via Josephson junctions 0| (see Figll3b). The central 
qubit (denoted by c in the figure) will encode the state 
to be cloned while the upper and lower qubits (denoted 



with u=up and ci=down) are initially in the blank state. 
All the Josephson junctions are assumed to be tunable 
by local magnetic fluxes. The total Hamiltonian of the 
3-qubit system is given by the sum of the Hamiltonians 
of the qubits Hq plus the interaction between them H cou . 



H 



]T 8E c af>-Ejaf 



(46) 



where Ej is the Josephson coupling in the Cooper pair 
box and 8E C is the energy difference between the two 
charge states of the computational Hilbert space. The 
coupling Hamiltonian for the 3-qubit system is 



i—u.d 

- (1/2) J2 + h - c -l (47) 

i— u,d 

Here Jk is the Josephson energy of the junctions which 
couple the different qubits and o± — [p x ±ia y )/2. If the 

coupling capacitance between the qubits is very small as 

(i) 

compared to the other capacitances one can assume E£ 
to be negligible. In practice, however the capacitive cou- 
pling is always present therefore it is necessary to have 
Jjc 3* 4Ejj . Then the dynamics of the system ap- 
proximates the ideal XY dynamics required to perform 
quantum cloning. The protocol to realize the SNC re- 
quires the preparation of the initial state. This can be 
achieved by tuning the gate voltages in such a way that 
the blank qubits are in |0) and the central qubit is in 
the state to be cloned. During the preparation the cou- 
pling between the qubits should be kept zero by piercing 
the corresponing SQUID loops of the junctsion Jk with 
a magnetic field equal to a flux quantum. In the second 
step, Hq is switched off and the dynamics of the system 
is entirely governed by H cou . At the optimal time the 
original state is cloned in the u and d qubits. 



0.854 




08 E,/J 



k k 



FIG. 14: Fidelity for the coupling Hamiltonian l|47|l as a 
function of Ek/Jk- 

As the implementation with superconducting nanocir- 
cuits has a slightly different Hamiltonian as compared to 
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the ideal XY model it is important to check for the loss 
of fidelity due to this difference. As it is shown in Fig |14l 
for Jk/Ek < 0.1 the maximum fidelity achievable differs 
at most by ~ 10~ 2 from the ideal value. 



X. CONCLUSIONS 

We have demonstrated that quantum cloning, in par- 
ticular PCC, can be realized using no external control 
but just with an appropriate design of the system Hamil- 
tonian. We considered the Heisenberg and XY coupling 
between the qubits and we found that the XY model 
saturates the optimal value for the fidelity of the 1 — ► 2 
PCC. In all other cases we have analyzed (N — > M PCC, 
universal cloning, cloning of qudits) our protocol gives a 
value of the fidelity of clones that is always within a few 
percent of the optimal value. As compared to the stan- 
dard protocol using quantum gates, however, there is a 
major advantage. Our setup is fast and, moreover, its ex- 
ecution time does not increase with the number of qubits 
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